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.:  Jig? 

In  this  thesis  the  tortional-flerural  response  of  a  t     mi* 
sional  airfoil  to  forced  oscillations  of  various  frequencies  and  at 
various  airspeeds  is  investigated.  The  airfoil  o'  oseci  has  character- 
istics vrhich  ore  typioal  of  modem  ^Tcrican  transport  wings a  and  the 
speeds  oover  the  range  fron  zero  airspeed  up  through  tho  speeds  for 
torcional-flexural  flutter  arid  torsional  divergence. 

In  the  latter  part  of  this  thesis,  curves  are  plotted  showing 
the  effects  v/hich  changes  in  the  assuned  wing;  parameters  have  on  tl\o 
torsianal-floxural  fluttor  speed  and  the  torsional  divor^orioe  speed. 


ctioi; 

study  of  the  phenomenon  of  flutter  of  an  airfoil  in  an  airstream 
leads  at  onoe  to  the  conclusion  that  it  is  very  complicated  in  nature 
and  consequently  practically  impossible,  or  at  least  very  difficult,  to 
visualize  clearly.  The  equations  of  motion  of  an  airfoil  osoi Hating  in 
an  airstream  have  been  derived  by  a  number  of  investigators;  and  the 
conditions  at  flutter,  or  at  the  limit  of  stability  of  the  wing,  have 
been  investigated.  In  this  thesis,  the  conditions  which  exist  at 
velocities  and  frequencies  other  than  flutter  are  investigated  in  hopes 
that  such  an  investigation  will  shed  additional  light  on  the  nature  of 
flutter  and  how  it  develops. 

The  type  of  flutter  selected  for  this  investigation  is  the  tor- 
si onal-flexural  type  or  one  which  involves  the  bending  and  twisting  of 
the  wing.  Such  a  system  is  one  of  two  degrees  of  freedom.  The  main 

orence  between  the  oscillating  wing  and  the  usual  type  of  oscillating 
system  that  we  are  accustomed  to  consider  is  that  the  wing  in  itself  is 
not  a  conservative  syctem,  for  under  the  proper  conditions  it  will  absorb 
energy  from  the  airstroora,  and  it  is  this  fact  that  makes  flutter  pos- 
sible. The  flexing  and  twisting  of  the  wing  in  the  airstream  cause  the 
aerodynamic  forces  and  momenta  on  the  wing  to  cluoige.  These  changes  in 
turn  cause  the  wing  to  flex  and  twist,  *n  general,  the  magnitude  and 
phase  relation  of  these  changes  are  functions  of  the  velocity  and  fre- 
quency. The  majority  of  the  time,  the  net  effect  of  these  forces  and 
moments  is  to  damp  the  oscillation,  but  under  oortain  conditions  their 
phase  relatione  may  bo  such  that  they  hare   a  negative  damping  effeot 


end  cause  the  oscillations  to  build  up.  This  condition  vve  call  flutter. 

Most  of  the  theoretical  investigations  of  flutter  have  beon  two 
dimensional,  that  is,  the  variations  of  the  parameters  along  the  vring 
span  have  been  neglected  or  at  the  best  averaged.  For  this  reason,  and 
because  the  parameters  involved  are  difficult  to  determine,  it  has  been 
suggested  that  it  may  be  possible  to  determine  the  flutter  speed  of  an 
actual  plane  by  plotting  ourves  of  the  natural  frequenoy  of  the  wing  in 
bending  and  torsion  against  airspeed  and  then  extrapolating  these  curves 
to  Intersection  under  the  supposition  that  the  speed  at  T.iiioh  the  two 
natural  frequencies  ooinoide  is  the  flutter  speed.  The  nature  of  the 
ourves  must  be  known  in  order  to  perform  this  e;:trapolation  and  will  be 
investigated  in  this  thesis*  The  natural  frequencies  could  be  determined 
in  flight  by  means  of  a  forcing  osoillator  and  pick-up  or  by  means  of  a 
pick-up  alone  allowing  the  accidental  "air  bumps"  to  set  the  wing  in  vi- 
bration. 


IIOTATIOIJ 


o  ■  cliord  of  wing 

C     ■  distance  of  the  olastio  axis  behind  t2ie  quarter  chord  point  in 

fraction  of  chord. 

a  ■  distance  of  the  centor  of  gravity  or  inertia  axis  behind  the  elas- 
tic axis  expressed  ac  a  fraction  of  the  chord. 

e  »  elastic  axis,  i.e.  the  point  at  whioh  an  applied  force  Trill  produce 

no  rotation  of  the  wing. 

k  -  spring  oonstant  of  the  wing  in  bending  of  flexure  (force/unit 

vertical  displacement), 

T/$*  torsional  constant  of  the  wing  (torque/radian). 

T/» 
^o  -   k   or  the  "elastic  radius" ^angular  dofleotion  of  the  wing. 

y*   •  vertical  displacement  of  any  point  j  of  the  wing,  positive  u;  , 

P  •  oircular  frequency  of  oscillation  (radians/sec). 


1/  x»f  natural  frequency  of  the  wing  in  bending  neglecting  the  effect 
of  the  apparent  raass  of  the  air. 


4. 


bl,  ■  nass  of  wing  per  unit  span. 

hit  ■  "apparent  mass  of  air",  i.e.  the  iaass  of  the  circumsoribed 

cylinder  of  air  one  unit  in  length. 

t  *  forward  velocity  or  airspeed. 

V  •  v  •  reduced  velocity. 
i>  o 

"p  •  ^-iB  ■  a  oomplex  function  of  the  reduced  velocity  "which  enters 

into  the  consideration  of  the  effect  of  the  wake 
on  a  wing  in  steady  state  oscillations. 

(  )  •  static  response,  tliat  is,  the  response  at  P   •  0  and  v  ■  0. 


pa 


Ph 

L 

I 


vertioal  osoillatory  force  applied  at  the  point  j,  positive  up. 

naraent  of  t}ie  applied  force  about  the  point  j,  positive  nosing  up. 

mass  density  of  air. 

.25  chord  point. 

.50  chord  point. 

^75  chord  point. 

lift  force 

aoraent  of  inertia. 


__  lUlIAL-FLEXURAL  0SCILIATIQIT3 

The  oace  to  be  considered  is  a  two-dimensional  one  where  the 
airfoil  is  assumed  to  have  uniform  charact or 1 sties  alone  the  span 
and  to  bo  infinitely  long. 

The  aerodynamic  forces  and  moments  on  an  airfoil  in  steady  state 
oscillation  are  as  follows! 

(l)  The  forces  and  moments  which  arise  from  the  fact  that  a  cer- 
tain amount  of  air  around  the  wing  also  must  bo  accelerated  when  the 
wing  is  accelerated.  Tho  "apparent"  mass  of  air  involved  when  the  wing 

is  accelerated  perpendicular  to  its  chord  is  equal  to  bit  •  ;i Per     or 

4 
the  circumscribed  cylinder  of  air.  The  "apparent"  moment  of  inertia 


of  air  which  is  the  amount  involved  in  angular  accelerations  of  the 

4 

T55~  ' 
(2;  The  changes  ir.  tho  lift  and  moment  due  to  vertical  velocity. 


wing  is  equal  to  It  -  n.oo 


velocity,  and  the  angular  rotation  of  the  wing  considering 
the  motion  as  quasi-steady,  that  is,  the  motion  is  considered  slow- 
enough  to  allow  the  circulation  or  lift  to  follow  the  changes. 

(3;  Tho  forces  and  mamoaits  which  arise  from  tho  effect  of  the 
rorticc         a  vring  in  steady  state  oscillation.  This  classification 
ic  t  ornan  and  ^ears  in  roferonce  (l;. 

considerations  result  in  the  following  forces  and 
.ts  on  an  airfoil  in  steady  state  osoillationc: 
At  25  per  co.it  of  the  chord 

L  -  R^CV2  (  +  - 


G. 


At  50  per  cent  of  the  chord 

L  n/0°2  (&  +  h) 

At  75  per   cent  of  tho  ohord 


L  ■  nPo^  t  $ 
4 

Moment  couple 

-"VST  • 

where  (  )  and  (**)  are  the  first  and  second  tine  derivatives.  These 
relations  have  been  derived  by  Lombard  in  reference  (2). 

In  addition  to  the  aerodynamic  forces,  the  following  elastic 
and  inertia  forces  Trill  acti 

(1)  The  moments  and  forces  arising  from  the  deflection  of  the 
elastic  axis  and  the  rotation  about  the  elaotio  axis. 

F  «  -ky#       and     11  •  -lc/^e^ 

(2)  The  foroes  and  moments  arising  from  the  inertia  and  moment 
of  inertia  of  the  wing's  mass. 

F  ^-sopysj      and     II  •  -I  4 
Using  these  foroes  and  moments,  Lombard  in  reference  (2)  has  set 
up  the  equations  of  equilibrium  for  a  two-dimensional  wing  in  steady 
state  oscillation.  The  derivation  assumes  no  structural  damping  and 
small  oscillations,  the  latter  condition  allowing  second  order  and 
higher  terms  to  be  neglected  in  certain  instances.  The  equations  are 
as  follows t 
(la)  ♦£  -  yh  ["♦  i^2-(l*l)  \?  Z+±l  ~i?r~\ 

+  Jo[+i\zll   -£)  ♦  (£♦  a-  1  -  l)?2-*?^   -U?xl 
L     2  2   fy         p   C"   \i     cj 


7. 


(lb)  ♦  M.250  -5hT-  *\2£*(£«-  a  +  _l)i>2[ 
amp       L  4u    J 


In  the  above  equations  the  barred  quantities  are  vectors  rotating 
in  the  eomplex  plane  at  an  angular  velocity  ^.  Tho  real  parts  of  these 
veotore  represent  the  actual  quantity  and  thoir  relative  angular  posi- 
tione  determine  their  pliase  relatione.  This  use  of  the  complex  plane  is 
cannon  practice  in  many  fields  of  engineering,  especially  in  electrical 
engineering,  and  it  is  particularly  adapted  to  vibration  problems.  A 
thorough  explanation  of  its  use  is  given  in  "Mechanical  Vibrations"  by 
J.  P.  Den  Hartog. 

To  determine  the  limit  of  stability  or  flutter  conditions  the  ap- 
plied foroes  and  moments  are  set  equal  to  zero. 

The  equations  then  have  the  formi 

0  -  yk  All  ♦  ♦«  A12 

0  "  Yh   A21  +  ♦«  A22 

In  order  that  there  be  a  solution  other  than  0,0  the  determinant 

of  the  coefficients  of  y  and  <J>o  must  equal  zero,  and  since  this  involves 
complex  quantities,  we  obtain  two  equations  to  determine  V  and  P from  the 
condition  tliat  both  the  reals  and  imasrinaries  must  equal  zero*  Since  r 

iot  be  expressed  explicitly,  the  equations  oannot  be  solved  explicitly, 
but  Ilasaner  and  Fingado  in  reference  (3)  have  provided  ourves  for  the  so- 
lution of  the  equations. 

This  solution  elves  only  the  conditions  at  flutter  speed  and  it  is 
desirable  to  investigate  the  conditions  which  exist  at  other  speeds  and 
frequencies. 


Because  of  the  number  of  parameters  whioh  enter  into  th#  equi- 
librium equations,  it  was  neoeGsary  to  assume  oertain  numerical  values 
for  all  of  them  exoept  the  velocity  of  flight  and  tlie   frequency  Of 
vibration.  In  choosing  these  numerical  values  an  attempt  was  made  to 
select  ones  which  were  typical  of  modern  American  transport  airplanes. 
The  values  selected  vrere  as  follows  i 

1«  Elastic  axis  at  35  per  cent  of  the  chord,  i.e.  £■  0.10 
2.  Center  of  gravity  of  v.lng  at  40  per  cent  of  chord,  i.e., 

£  ♦  c  •  G.15 


3.  Radius  of  gyration  of  wing  «y0. 0625c2  i.e.  if  •  0.0625 

4.  Bonding  frequency  at  zero  airspeed,  neglecting  the  apparent 
mass  of  tile  air 

»  300  eyeles/nin  *  10n  rad/seo 
5;  elastic  radius  »^-.5C  i.e.,   AJ2  -  0.50 

6»  Ratio  of  the  mass  of  wing  per  unit  length  to  the  apparent  mass 
of  the  air  H  »  6.0C 

7.  fling  chord  «  7.5  ft. 

8.  Y/ing  weight  -  2.7  lbs/ft 

Then  by  Kacsner  lingado  Chart  flutter  oocurs  at$ 
&       •  57.2  rad/sec 
v/o   »  73.2  l/soo  corresponding  to  374  mph. 
V   *  1.28 
burning  now  to  the  equations  of  equilibrium  (la)  h   (lb)  and 
substituting  the  wing  parameters,  we  obtain i 


9. 

An  s  , A,- 


(2a)  +  7_  -  y.   [+  987  ♦  0.667  i^r  7-1.167^  2J   +  *?o  [+395-0.667  t2  7 

-0.167  i^v-  0.3917P»2J 

, Aii ,  , .  Azt , 

(2b)   +  l!.25o-  yh  ["-98.7  ♦  0.1917  P^]   ♦  ♦©  [♦  454-0.OO48^2+  0.0633  i  P r  1 

-.otei     Thic  notation  of  A,,#  A^g,  etc.  does  not  agree  exactly  with  the 
notation  of  -^i\$  &\2»  e^c#  ^n  kombard'8  report. 


Dsr^riiMATiou  on'.T^,  : .: ■■:;-\l  ltod::s  at  t»o 

Setting  v  *  0  in  equations  (2)  we  obtains 
(3a)     +  7_  -  y^  [+  907  -  1.1G7  P2]+  To  [♦  396  -  0.3917 PZJ 

(3b)     ♦  H.25o  -  y^T-  98.7  ♦  0.1917P2]*-  $e  [  +  454  -  0.0048^ 

c:  >;, 

To  determine  the  natural  frequencies  of  vibration  we  set  the 
determinant  of  the  coefficients  of  y  and  $0  to  zero  and  obtain  the 
frequency  equations 

p4  -  8030  0  2  ♦  6,030,000  -  0 
the  roots  of  ..v.ich  ares 
P   *  29.0  rad/sec 
P   -  34.8  rad/seo 
Nov  using  equation  (3b)  and  substituting  K»  29.0,  v*s  obtain 

0  *  Jh  C"  9-7   1C1]  *  *°  [*  454  -  *»°] 
<>e  •  -.138  y^ 

This  equation  in  satisfied  if  the  wing  oscillates  with  a  nodal  point 
6.5  sherd  lengths  ahead  of  the  loading  odre. 


10. 

Again  using  equation  (ob)  and  substituting  P  »  84.8,  \ro   obtain 
0  •  yh  [-98.7  +  1378]  +  *?o  [+  454-34.5] 

*$c  -  -3.06  yh 

This  equation  indicates  a  nodal  point  at  42.2  per  cent  of  the  chord 
behind  the  leading  edge. 

If  an  oscillating  force  ic  applied  at  tho  nodal  point  of  one 
normal  node,  the  resulting  vibrations  will  be  in  the  other  nornal  aoda, 

Kt  v/ill  nor/  shift  the  flexural  coordinate  of  our  equations  from 
yl  to  y  422o*  •^ien  froia  the  geometry  of  the  set-up,  v/e  obtain 

yk  •  7;422  -   •528  "^ 
cm  substituted  in  equations  (2)  results  int 

(4b)  +  A^  .  ^.,, ,  p^  J  +$c[A^-.3k*  A\  J 

1'urther  transformation  will  plaoe  these  equations  in  a  nore 

useful  form.     Noting  that 

ft._     -  H  or     -7x  G.75c 
-C.5o  -o 


and  substituting  (4a)  and  (4b)  in  this,  ire   obtaint 

r = (?•*■  *  /^/  ~ *• ^J + ^7  4 ;. -  ^ ? «£,  -cists z **$} 


c  m  & 
or 


Likewise  j 

B  _0_     .  S  9-     +Fx  .172c     which  when  (4a)  and  (4b)  are 
substituted  in,  results  in 


(5b 3    ' ^^  -  ^«  //£ ,  +. 7 ^  ^ ]*4 c[-  c i,L  +A„  t.PZ  Ar<  - .  32  *'4 /  4  j 


11. 


oting  that 

G.75  m*,*       6#75cjap 

M1^  g*^~*  -l.oc  -  ",422c  f  equations  (5)  then  reduce  to  the  following 
Dp        amp 

■whan  the  values  of  the  A»s  from  (2)  ore  Introduced, 


(6a)     L°**E-  +  i*c  =  #«*<.  [+  tOOZ  +0.<cG7+PPvyc  -  /./75"P*7 


*  <pc  V-0,z{  tf^i  pp*r/z-  cm < v4/c  -  O.G>t>7P«ryc*J 


a 


_  ^- ■ "zi ^ 

(6bj     6.9£  /--d.^c  wji^S+V.I  +  O.I t+7+  *> P*r/c  -.  0090^] 


B: 


~>2Z— > 

Equations  (C)  are  not  actually  In  the  normal  form  at  sero  velocity, 

because  the  normal  coordinates  are  y  i0<5.  and  y  e  r   ,  but  since  the  lat- 

•/.4<;2o     j-6.5q' 

tor  point  does  not  lie  on  the  wing,  and  since  it  is  desirable  to  retain  $ 
as  a  coordinate  in  order  to  better  visualise  the  notion,  the  equations 
have  boon  placed  in  this  form. 

The  equations  of  motion  are  now  such  that  we  can  obtain  the  two 


principal  nodes  of  oscillation  at  zero  airspeed  by  setting  either  F  A99 

or  F  ..  equrd  to  zero,  thus  applying  the  oscillating  force  at  a  nodal 
— o#oc 

point,  The  result ing  oscillations  will  be  about  the  opposite  nodal  point. 


i:odu  i 

To  produoe  what  we  sliall  call  Mode  1,  we  apply  an  oscillating 
forco  of  the  form  F'coaPt,  or  in  the  oomplex  plane  of  the  form  r'e   , 
at  the  ,422  chord  point,  and  in  order  to  express  the  amplitude  of  re- 


;i2 


aponse  In  a  form  whioh  is  independent  of  the  magnitude  of  the  applied 
force,  wo  divide  by  the  static  deflection  which  this  applied  force 
would  produoe  if  the  airspeed  (▼)  and  the  frequency  (*>)  were  tero. 
Equation*  (6)  then  become 


(7b)    6.9ZA+sSm     Q     =g.„z[B2,]  +  tc[B2J 


(7o)     £7, 


-9ZZ        +  — r  € 


2L 


.42Z 


£>./a   e-""Jt 


(9-<~K     (&*M) 


(7d)      <p 


B. 


c  =    —  -=■ 


// 


z.\ 


Fe 


*>t 


(7«)    ^e 


&/ 


^« 


wher»    ^ 


and  (   )#  denotes  the  static  deflection. 
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HOPE  Z 
In  a  eiiailar  nanner  we  can  produce  Mode  2  by  applying  the  osoil« 
latin£  force  at  the  -6.5  chord  point, 

(8a)     /.OZ5F.^t  m        Q        __   g^zz  £  £J    +  4c[B,k  ] 

(3b)    £>^Zj-L.s*  .  F>e;*tm    -^  j-  q^ j  +  ^  j-  g£J 


(3o) 


(8d) 


-J.*  '^Z 


<fa    =     + 


B*-<.  p\,-ot 


z\ 


fe 


A/  f'^** 


A 


-i«x;i 


^A     £>*>* 


#o.     fe,/^x 


&<0>  "    (B,.//x)o 


BJa  e*** 


(8.)    3>»**   .  _  _&* 

where  A  and  (   )#  have  the  some  significance  as  in  Mode  1, 

In  equation  (8o)f  the  flexural  response  has  boon  divided  by  the 

static  angular  response  because  the  static  flexural  response  of  the 

.422o  point  is  zero. 


:::  or  :  -.;  -t.'g  at  various  velocities  Ain) 

FRK^lTJIICISS 
e  above  equations  enable  us  to  obtain  response  ourves  for  the 
ig  at  any  velocity,  Tho  procedure  followed  in  obtaining  the  response 
curves  in  thi3  thesis  -rcas  firct  to  choose  a  value  of  v/o.  This  was 


14, 


substituted  in  equations  (7)  and  (8).  Values  of  p  from  0  to  120 
in  steps  of  10  v/ere   then  assumed  and  the  value  Y  determined  from 
the  corresponding  V  «  v/pc.  This  enabled  us  to  caloulate  ths  B's 
and  henoe  the  deflections.  The  calculations  are,  of  course,  in 
the  complex  form.  The  deflections  plotted  are  the  maximua  deflec- 
tions j  that  is,  the  modulus  of  y  and  ^. 


15, 

TORSIONAL  DIVl2taEHCE 

■■■      ■»     ■  —■■■■■    ^m      m  ■■     ■     m.  m    —  ».  ■      ■     ■■  m 

Torsional  divergence  will  ooour  in  the  wing  when  the  rata  of 
ohange  of  moment  of  the  aerodynamic  forces  with  re.  L>eot  to  a  change 
in  the  angle  of  attack  exceeds  the  rate  of  ohange  of  the  elastic  re- 
storing moment  with  respect  to  a  ohange  in  angle* 

In  order  to  determine  the  velocity  at  whioh  this  divergence 
occurs,  we  must  deternino  the  velocity  at  whioh  A  of  our  equations  of 
motion  bee  ones  equal  to  sero  at  zero  frequency,  since  this  is  the  fre- 
quency at  which  divergence  occurs* 

Then  noting  that  ati^«  0, 


becomes 


V  •  v_  moo  end  V  *  \9  the  A  of  equation  (la)  and  (lb) 
t>  e 


[P,*J        t>f('/z-£)-±*g] 


fP/<J      [-*>,*  {ei&-e)-Z 


z) 


'_J 


*  o 


from  r/hioh  we  obtain 


rd  fa  •    J^JI  iLjL.      Rt  divergence 


e 


Substituting  our  parameters 

Tj  fa  ■  1°*  x  ,707    /  6     -  86  rad/sec 

d        — g y  j 


v   •  86  x  7.5  -  646  ft/sec  -  440  a.p^h. 


Response  curves  for  this  velocity  have  been  plotted. 
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.,.;  .••_.>:,  £__  ;ns  ..i:;c-  . : j.  .::.  _^___  0£  r;6.;i;:xs.i 

©  can  obtain  on  equation  for  the  total  applied  force  on  the 
wing  by  adding  the  foroo  at  .422o  to  the  force  at  -6.5o.  Thon  by 
using  equations  (6),  vd  obtain 

<rtp  Tip  *T7fi 

^ C, .  ^- C>  i:  — -~> 

<9a'     b  -  -    -?  r    a  S" 


a. 


If  ire  now  set  "Jo  equal  to  zero,  to  will  obtain  an  equation  for 
tl»  motion  of  a  vring  with  infinite  torsional  rigidity.  By  setting 
$o  «•  0,  we  imply  tiiat  we  are  applying  the  force  so  tliat  its  moment 
varies  in  such  a  maimer  tliat  the  moment  equation  (Ga)  or  (Gb)  is  also 
satisfied  for  $e  ■  0. 

.uation  (9a)  is  equivalent  to  equation  (4a),  but  as  tli©  B'a 
have  already  boon  calculated,  it  is  the  more  convenient  form, 

-•ow  again  applying  an  oscillating  force  of  the  form  F'  oos^t 

iPt 
or  F'e  '   in  the  oomplex  plane  and  dividing  the  response  amplitude 

by  the  static  response,  wa  obtain 

<»)    _£~Iit  -  F'eT'*-  #,«*/'£/ 7  *■  o[c/z] 


/7>f. 

(9o) 

§•*"■ = 

r-'       Atft 

r  e 
£44 

(M) 

^g  ,+*  x. 

(ctt)o 

jS«***Jo  ^-//  (flexure  only) 


17, 


Uoturning  now  to  the  equation  (6b)  and  setting  y     equal  to 
:;ero,  we  obtain  the  equation  of  notion  of  the  wing  with  the  .422  chord 
'■joint  fixed. 


(10a)      '?£  f- 


o[BZl]  +$c[B„] 


It  is  alco  implied  that  we  are  applying  a  foree  at  the  .422  chord 
point  of  such  variable  magnitude  that  equation  (6a)  is  satisfied  for 
the  condition  y  ^2  m  ®* 

Then  applying  an  oscillating  inoment  of  the  form  H'cos  Pt,  or 

i  ^t 
!!••     in  the  complex  plane,  and  dividing  the  response  by  the  statio 

response,  we  obtaini 


(10b)     <£c  ■ 

A4'      <'t>t 

A1  e 

5u 

(lOo)     <#C 

(8  a  4.)  O            -<%*>t 

-    — _ e 

(&C)  Pi  (*422  chord  point  fixed) 

:  zponse   curves  for  the  above  two  modes  of  oscillations  have 
been  plotted.  Tie  amplitudes  plotted  are  the  maximum  values  or 
modulus  • 
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DISCUSSIOII  OF  THE  BSSPOKSE  CURVES 

The  response  curves  for  one  degree  of  freedom  as  shmvn  in 
Figures  1  and  2  are  typieol  families  of  damped  oscillations  -with  the 
damping  increasing  with  airspeed.  The  torsional  response  curves  show 
the  effeot  of  the  torsional  divergence  phenomenon  near  zero  frequoncy 
with  the  response  boooming  large  at  332  m.y.h.  which  is  slightly  below 
the  divergence  speed  of  the  wing  with  the  .422  chord  point  fixed.  At 
speeds  above  this,  the  deflection  is  actually  opposite  to  the  applied 
force  which  means  that  a  foroe  must  bo  applied  to  keep  tlie  wing  from 
diverging. 

The  response  curves  given  in  Figures  3  to  11  for  **ode  1,  or  the 
node  which  at  zero  airspeed  is  predominately  flexural,  indioate  that 
natural  frequency  of  flexure  which  occurs  at  zero  airspeed  is  damped 
out  as  the  airspeed  increases.  The  damping  is  even  greater  in  the 
two  degrees  of  freedom  than  in  the  one  degree  of  freedom  case  discus- 
sed above.  It  should  be  remembered  in  examining  these  curves  that  as 
soon  as  the  airspeed  is  other  than  zero,  the  modes  are  no  longer  nor- 
mal modes,  and  we  are  obtaining  some  cross  responses.  The  peaks  on 
the  right  of  the  curves  of  response  of  y  .„  to  a  foroe  at  the  .422 
ohord  point  are  of  this  nature  since  they  arise  from  the  $o  or  tor- 
sional response. 

The  curves  of  Mode  2  or  the  torsional  mode  at  zero  airspeed,  as 
given  in  Figures  12  to  19,  indioate  tliat  the  torsional  oscillations  are 
increasingly  damped  with  increasing  airspeed  and  that  the  frequency  of 
maximum  response  movos  to  the  left  or  decreases.  However,  as  the 
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flutter  speed  i8  approached,  the  damping  bo  com©  a  less  and  the  amplitude 
builds  up.  The  curves  for  440  m.p.h.  indicate  torsional  divergence  oc- 
cur in,-  at  zero  frequency. 

As  stated  in  the  introduction,  one  of  the  purposes  of  this  inves- 
tigation was  to  determine  the  possibility  of  predicting  the  speed  at  which 
flutter  would  oocur  by  plotting  the  natural  torsional  and  flexural 
frequencies  at  various  airspeeds  as  determined  by  flight  tests.  If  the 
frequency  at  which  the  maximum  res pons©  occurs  at  any  given  airspeed  is 
called  the  natural  frequency  and  the3e  frequencies  plotted  against  air- 
speed, the  ourves  will  be  those  whioh  would  be  obtained  by  flight  tests 
and  are  shonn  in  Pi£ure  21.  It  can  be  seen  at  once  that  flutter  does 
not  occur  where  the  torsion  and  flexural  curves  intersect,  but  that  the 
flexural  mode  has  been  damped  out,  as  was  pointed  out  before,  and  that 
flutter  occurs  at  a  point  on  the  torsional  curve  without  warning  other 
than  a  building  up  of  the  amplitude  of  response.  This  method  of  attempt- 
ing to  predict  flutter  appears  to  be  irapraotioal. 


FIGURE   1 
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FIGURE  3 
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effect  of  ciiakoss  r,  tie  various  \aiio  parameters 

In  order  to  obtain  soae  idea  of  the  importance  of  the  various 
parameters  and  to  investigate  the  effects  which  changes  in  these  para- 
meters -.rill  have  on  the  torsional  flexural  flutter  speed  and  the  tor- 
sional divergence  speed,  one  parameter  at  a  tine  ha3  been  varied  while 
holding  the  remaining  wing  parameters  at  the  values  chosen  for  the  re- 
sponse investigation,  "y  use  of  the  Kaesner  and  Fingado  chart  of 
reference  (3)  for  flutter  speed  and  the  previously  derived  equation  for 
torsional  divergence  speed,  curves  have  been  plotted  showing  thoir 
variation  and  are  given  in  Figures  22  to  29,  The  ciroled  points  in- 
dicate  the  basic  wing. 

Tlxe  torsional  divergence  being  non-dynamic  in  nature  is  rather 
easy  to  visualise.  Obviously  the  bending  frequency  and  inertia  axis 
have  nothing  to  do  with  torsional  divergence  henoe  the  torsional  di- 
vergence speed  remains  constant  as  they  are  varied.   Increasing  the 
wing  weight  while  holding  the  torsional  frequency  constant  or  increas- 
ing the  torsional  frequency  while  holding  i  and  bl,,  i.e,  tho  moment  of 
inertia,  constant  means  that  the  torsional  stiffness  and  hence  the 
divergence  speed  is  increasing.  Moving  the  el actio  axis  back  on  the 
wing  in  effect  gives  the  aerodynamic  forces,  which  oan  be  considered 
as  acting  at  the  quarter  cliord,  a  longer  moment  arm  and  henoe  decreased 
the  divergence  spoed.  Of  course,  the  elastic  axis  at  the  quarter  chord 
point  would  result  in  an  infinite  divergence  speed.  An  increase  in 
altitude  decreased  the  aerodynanio  forces  on  the  wing  at  a  given  speed 
and  thus  increases  the  divergence  speed. 
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The  reasons  bcliind  the  changes  in  the  flutter  speeds  are  very 
complex  and  on  explanation  will  not  be  attempted,  *-»owever,  the  follow- 
ing pertinent  points  should  be  noted s 

1.  The  flexural  frequenoy  or  stiffness  is  apparently  of  minor 

inp  orb  once  • 

2.  The  flutter  speed  increases  vrlth  increased  torsional  frequenoy 

in  nearly  a  linear  relation. 

3.  Moving  the  elastic  axis  to  the  rear  increases  tho  flutter  speed. 

4.  The  position  of  the  inertia  axis  has  a  marked  effect  on  the 

flutter  speod,  a  forward  position  being  favorable. 

5^  From  the  curve  showing  the  effect  of  changes  in  the  radius  of 
gyration  of  wing,  it  should  be  noted  that  even  at  zero 
radius  which  corresponds  to  a  high  torsional  frequenoy  the 
flutter  speed  is  still  quite  nominal.  From  this  we  can 
conclude  that  the  major  effect  of  inoreasod  torsional  fre- 
quency comes  from  an  inerease  in  torsional  stiffness. 

6^  Changes  in  wing  weight  have  a  peculiar  effect.  Since  the  fre- 
quencies have  been  hold  constant  part  of  the  effect  can  be 
attributed  to  the  fact  that  the  stiffnesses  must  change  as 
the  "weight  is  varied. 

7.  The  flutter  speed  increases  with  altitude. 

8.  It  should  bo  remembered  that  a  definite  wing  has  been  oonsiderod 

and  that  only  one  parameter  has  been  varied  at  a  time.  If 
another  wing  were  considered,  or  if  two  or  more  parameters 
were  varied  at  one  time  the  changes  in  flutter  speed  and 
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even  their  character  might  -well  be  different.  It  will 
be  noted  that  the  effect  of  moving  the  elastic  and  inertia 
axes  simultaneously  is  not  the  sum  of  the  effects  of  moving 
than  individually. 
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